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>— «^ ■ It is shown that most of the models for analyzing meson-baryon reactions in the 

^*_i ' micleon resonance region can be derived from a Hamiltonian formulation of the 

^-^ ' problem. An extension of the coupled-channel approach to include nirN channel 

is briefly described and some preliminary results for the Af*(1535) excitation are 
n»«J ■ presented. 
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1. Introduction 

With very successful experimental efforts in the past few years, we are now 
facing a challenge to interpret very extensive data of electromagnetic meson 

P^ . production reactions in terms of the structure of nucleon resonances (N*). 

To achieve this goal, we need to perform amplitude analyses of the data in 
order to extract N* parameters. We also need to develop reaction models 
•Oj . to analyze the dynamical content of the extracted N* parameters. At the 

C$ ' present time, we can use the N* data to test the predictions from various 

QCD-based hadron models. In the near future, we hope to understand N* 
parameters from Lattice QCD. 

In the A region, both the amplitude analyses and dynamical reaction 
models have been well developed. We find that these two efforts are com- 
plementary. For example, the 7iV — ► A Ml transition strength extracted 
from all amplitude analyses is Gm(0) = 3.18 ± 0.04 which is about 40 
% larger than the constituent quark model prediction. This difference is 
understood 1 ' 2 by developing dynamical reaction models within which one 
can show that the discrepancy is due to the pion cloud which is not included 
in the commonly considered constituent quark model prediction. 

In the second and third resonance regions, the situation is much more 
complicated because of many open channels. It is necessary to develop 
coupled-channel approaches for learning about the N* properties. The 



main objective of this contribution is to review the development in this 
direction. We will also describe a newly developed coupled-channel model 
which is aimed at accounting for rigorously the ttttN unitarity condition. 

In section 2, we will introduce a Hamiltonian formulation within which 
most of the current models of electromagnetic meson production reactions 
can be derived and compared. The extension of the coupled-channel ap- 
proach to account for explicitly the ttttN channel is then described in section 
3. A summary is given in section 4. 

2. Derivation of Models 

The starting point of our derivation is to assume that the meson-baryon 
(MB) reactions can be described by a Hamiltonian of the following form 

H = H a + V, (1) 

where Ho is the free Hamiltonian and 

V = v b9 + v R . (2) 

Here v bg is the non- resonant (background) term due to the mechanisms such 
as the tree-diagram mechanisms illustrated in Fig. l(a)-(d), and v R de- 
scribes the N* excitation Fig. 1(e). Schematically, the resonant term can 
be written as 

N* l 

where Ti defines the decay of the i-th N* state into meson-baryon states, 
and M® is a mass parameter related to the resonance position. 

The next step is to define a channel space spanned by the considered 
meson-baryon (MB) channels: jN, irN, rjN, ttA, pN aN, ■■. The S-matrix 
of the meson-baryon reaction is defined by 

S(E) a>b = 5 a . b - 2mS(E - H )T a>b (E) , (4) 

where (a, b) denote MB channels, and the scattering T-matrix is defined 
by the following coupled- channel equation 

T a , b (E) = V afi + J2 V a .c9c(E)T c . b (E) . (5) 

C 

Here the meson-baryon propagator of channel c is 
g c (E) =< c | g(E) \ c > 
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Figure 1. Tree diagrams for meson- baryon reactions. N* is a nucleon resonance state, 
with 

9 ™ - E-H + i e 

= g p (E) - inS(E - H ) , (6) 



where 



9 P (E) = 



P 



Here P denotes taking the principal-value part of any integration over the 
propagator. If g{E) in Eq.(5) is replaced by g p (E), we then define the 
K-matrix which is related to T-matrix by 



T a . b (E) = K a . b {E) -J2Ta,c(E)[nrS(E - H )] c K c , b (E) 



(8) 



By using the two potential formulation, one can cast Eq.(5) into 

T a , b (E)=t% b (E)+% b (E) (9) 

with 

*&(£)= E r\ f . ta (E)[G(E)] itj r N . tb {E). (io) 



N'N't 



The first term of Eq.(9) is determined only by the non- resonant interaction 
t\(E)=vl% + Y J < 9 ME)t b c %{E). (11) 



The resonant amplitude Eq.(lO) is determined by the dressed vertex 

T N . >a (E)=T N . ia + J2 T N;b9b(E)tb! a (E), (12) 

b 

and the dressed propagator 

[G(E)-\ 3 (E) = (E- M° N: )S h] - Z id (E) . (13) 

Here Mjy, is the bare mass of the resonance state TV* , and the self-energy 

is 

ZiAE) =J2 T N*,a9a( E ^N;AE) . (14) 

a 

Note that the meson-baryon propagator g a (E) for channels including 
an unstable particle, such as 7rA, pN and aN ', must be modified to include 
a width. In the Hamiltonian formulation, this amounts to the following 
replacement 

94E)^<a\ E _ Ho 1 _ Jlv{E) \a>, (15) 

where the energy shift is 

5V(£) = J2 r+(») P ™ N . Y v {i) . (16) 

i 

Here IV describes the decay of p, a or A in the quasi-particle channels. 

Eq.(5), Eqs.(9)-(16), and Eq.(8) are the starting points of our deriva- 
tions. From now on, we consider the formulation in the partial- wave repre- 
sentation. The channel labels, (a^bjc), will also include the usual angular 
momentum and isospin quantum numbers. 

2.1. Unitary Isobar Model (VIM) 

2.1.1. MAID 

The Unitary Isobar Model developed 3 by the Mainz group is based on the 
on-shell relation Eq.(8). By including only one hadron channel, irN (or r/N 
), Eq.(8) leads to 

TkN^n = e l5 ' N cosS 7:N K 7TNnN . (17) 

where S^n is the pion-nucleon scattering phase shift. By further assuming 
that K = V = v b9 + v R , one can cast the above equation into the following 
form 

T nN , 7N ( UIM) = e s ^cosS wN [v% nN ] + ]T T$„ N (E) . (18) 

N* 



The non-resonant term v bg in Eq.(18) is calculated from the standard 
Born terms but with an energy-dependent mixture of pseudo-vector (PV) 
and pseudo-scalar (PS) ttNN coupling and the p and u> exchanges. For reso- 
nant terms in Eq.(18), MAID uses the following Walker's parameterization 5 

t ^n^n( e ) = Hn{ e ) m? „ ° E 2 1 iMiT 101 ^jv(^)^' ' ( 19 ) 

where fl N {E) and f* N (E) are the form factors describing the decays of 
N* , Ttot is the total decay width, A 1 is the 7^ — ► N* excitation strength. 
The phase $ is required by the unitary condition and is determined by an 
assumption relating the phase of the total production amplitude to the ttN 
phase shift and inelasticity. 



2.1.2. JLab/Yeveran UIM 

The Jlab/ Yerevan UIM 4 is similar to MAID. But it implements the Rcggc 
parameterization in calculating the amplitudes at high energies. It also uses 
a different procedure to unitarize the amplitudes. 



Both MAID and JLab/Yeveran UIM have been applied extensively to 
analyze the data of 7r and r\ production reactions. Very useful new infor- 
mation on N* have been extracted. 



2.2. Multi-channel K-matrix models 

2.2.1. SAID 

The model employed in SAID 6 is based on the on-shell relation Eq. (8) with 
three channels: jN, irN, and irA which represents all other open channels. 
The solution of the resulting 3x3 matrix equation can be written as 

T 7 n,-kn(SAID) = Ai(l + iTjrJV.Triv) + AbT^n^n , (20) 

where 

a T s KjN,kaKkN,ttN / 01 n 

Ai = K lN ^ N — , (Zi.) 

-K-nN,TrA 

A R = f^ . (22) 



In actual analysis, they simply parameterize Ai and An as 

M 



Ai = [t^ )7rJV ] + Yl PnzQi a+n (z) , (23) 



n=0 

N 



K KO ^— J, "V 



where fco and go are the on-shell momenta for pion and photon respec- 
tively, z = \/k.Q + Am^/ko, Ql(z) is the legendre polynomial of second 
kind, E v — E~ ( — m w (l + m 7r /{2mN)), and p n and p n are free parameters. 
SAID calculates v 9 N ^ N of Eq.(23) from the standard PS Born term and p 
and u> exchanges. The empirical ttN amplitude T^n^n needed to evaluate 
Eq.(20) is also available in SAID. 

Once the parameters p n andp„ in Eqs.(23)-(24) are determined, the N* 
parameters are then extracted by fitting the resulting amplitude T-^^n 
at energies near the resonance position to a Breit-Wigner parameteriza- 
tion(similar to Eq.(19)). Very extensive data of pion photoproduction have 
been analyzed by SAID. The extension of SAID to also analyze pion elec- 
troproduction data is being pursued. 

2.2.2. Giessen Model 

The coupled-channel model developed by the Giessen group 7 can be ob- 
tained from Eq.(8) by taking the approximation K = V. This leads to a 
matrix equation involving only the on-shell matrix elements of V 

T a , b (Giessen) = £}(1 + iV{E))-\ tC V cjb {E) . (25) 

c 

The interaction V = v bg + v R is evaluated from tree-diagrams of various 
effective lagrangians. The form factors, coupling constants, and resonance 
parameters are adjusted to fit both the 7riV and 7^ reaction data. They 
include up to 5 channels in some fits, and have identified several new N* 
states. But further confirmations are needed to establish their findings 
conclusively. 

2.2.3. KSU Model 

The Kent State University (KSU) model 8 can be derived by noting that the 
non-resonant amplitude t bg , defined by a hermitian v hg in Eq.(ll), define 



a S-matrix with the following properties 

S b a 9 b {E) = S a , b 2mS(E H,)t%(E) (26) 

= J2^ C )T (E)^(E), (27) 

C 

where the non-resonant scattering operator is 

^( E ) = <c + 9a(E)t b a %(E). (28) 

With some derivations, the S-matrix Eq.(4) and the scattering T-matrix 
defined by Eqs.(9)-(14) can then be cast into following form 

S a AE) =J2^ T (E)RcAE)^(E) , (29) 

c,d 

with 

R c ,d{E) = 6 c , d + 2iT R d {E). (30) 

(31) 

Here we have defined 

T c R d (E) = ^2r^ <c (E)[G(E) iJ T Nhd (E) • (32) 

The above set of equations is identical to that used in the KSU model of 
Ref. 8 . In practice, the KSU model fits the data by parameterizing T R as 
a Breit-Wigner resonant form T R ~ xT/2/(E — M — iT/2) and setting 
c<j( + ' = B = B1B2 ■ ■ ■ B n , where Bi = exp(iX Aj) is a unitary matrix. 

The KSU model has been applied to nN reactions, including pion pho- 
toproduction. It is now being extended to investigate KN reactions. 

2.3. The CMB Model 

A unitary multi-channel isobar model with analyticity was developed 9 in 
1970's by the Carnegie-Mellon Berkeley(CMB) collaboration for analyzing 
the irN data. The CMB model can be derived by assuming that the non- 
resonant potential v hg is also of the separable form of v R of Eq.(3) 



v \,a T L,b r Ha T Htb 



(33) 



a - b E-M L E- M H 

The resulting coupled-channel equations are identical to Eqs.(9)-(16), ex- 
cept that i 9 6 = and the sum over TV* is now extended to include these 
two distance poles L and H. 



By changing the integration variables and adding a substraction term, 
Eq.(14) can lead to CMB's dispersion relations 

C 

Re[Ms)} = Re[Mso)} + °—^ f° M*c(*')] 

77 Js th {s' -s){s ' -so) 

Thus CMB model is analytic in structure which marks its difference with 
all K-matrix models described above. 

The CMB model has been revived in recent years by the Zagreb group 10 
and a Pittsburgh- ANL collaboration 11 to extract the N* parameters from 
fitting the recent empirical ttN and jN reaction amplitudes. The result- 
ing N* parameters have very significant differences with what are listed by 
PDG in some partial waves. In particular, several important issues concern- 
ing the extraction of the N* parameters in Sn channel have been analyzed 
in detail. 



2.4. Dynamical Models 

A. In the A region 

Keeping only one resonance N* — A and two channels a,b — irN, jN, 
Eqs.(9)-(14) are reduced to what were developed in the Sato-Lee (SL) 
model 1 . In solving exactly Eqs.(9)-(14), the non-resonant interactions 
v itN ttN an d v tt 9 n t-jv are derived from the standard PV Born terms and 
p and u> exchanges by using an unitary transformation method. 

In the Dubna-Mainz-Taiwan (DMT) model 2 , they depart from the for- 
mulation Eqs.(9)-(14) by using the Walker's parameterization defined by 
Eq.(19) to describe the resonant term t R of Eq.(9). Accordingly, their def- 
inition of the non-resonant amplitude also differs from Eq.(ll): t c 9 b in the 
right-hand side of Eq.(ll) is replaced by the full amplitude T c j,. Further- 
more, they follow MAID to calculate the non-resonant interaction vJ N N 
from an energy-dependent mixture of PS and PV Born terms. 

Extensive data of pion photoproduction and electroproduction in the A 
region can be described by both the SL and DMT models. However, the 
extracted 7^ — > A form factors, in particular their bare form factors, are 
significantly different. 

B. In the second and third resonance regions 

Eqs.(9)-(16) are used in a 2-N* and 3-channel (irN, rjN, and it A) 



study 12 of t:N scattering in Sn partial wave, aiming at investigating 
how the quark-quark interaction in the constituent quark model can be 
determined directly by using the reaction data. Eqs.(9)-(16) are also 
the basis of examining the N* effects 13 and one-loop coupled- channel 
effects 14 on lo meson photoproduction and the coupled-channel effects on 
K photoproduction 15 . 

The coupled-channel study of both ttN scattering and 7^V — ► ttN in 
Sn channel by Chen et al 16 includes irN, r)N, and 7^ channels. Their 
ttN scattering calculation is performed by using Eq.(5), which is of course 
equivalent to Eqs.(9)-(14). In their -yN — > irN calculation, they neglect 
the 7iV — > i]N — > ttN coupled-channel effect, and follow the procedure of 
the DMT model to evaluate the resonant term in terms of the Walker's 
parameterization (Eq.(19)). They find that four N* are needed to fit the 
empirical amplitudes in Sn channel up to W = 2 GeV. 

A coupled-channel calculation based on Eq.(5) has been carried out by 
Julich group 17 for ttN scattering. They are able to describe the irN phase 
shifts up to W — 1 .9 GeV by including irN, rjN, irA, pN and aN channels 
and 5 N* resonances : P 33 (1232), 5n(1535), Sn(1530), 5n(1650) and 
Di3(1520). They find that the Roper resonance Pn(1440) is completely 
due to the meson-exchange coupled-channel effects. 

A coupled channel calculation based on Eq.(5) for both irN scatter- 
ing and 77V — ► ttN up to W — 1.5 GeV has been reported by Fuda and 
Alarbi 18 . They include 7riV, 7^, r/N, and irA channels and 4 N* reso- 
nances : P 33 (1232), Pn(1440), 5n(1535), and Di 3 (1520). The parameters 
are adjusted to fit the empirical multipolc amplitudes in a few low partial 
waves. 

Much simpler coupled- channel calculations have been performed by us- 
ing separable interactions. In the model of Gross and Surya 19 , such sepa- 
rable interactions are from simplifying the meson-exchange mechanisms in 
Figs l.(a)-c) as a contact term like Fig. 1(d). They include only wN and 
7^V channels and 3 resonances: P 33 (1232), Pn(1440) and £>i 3 (1520), and 
restrict their investigation up to W < 1.5 GeV. To account for the inelas- 
ticities in Pn and -Di 3 , the N* — » irA coupling is introduced in these two 
partial waves. The inelasticities in other partial waves are neglected. 

A similar separable simplification is also used in the chiral coupled- 
channel models 20 ' 21 for strange particle production. There the separable 
interactions are directly determined from the leading contact terms of SU(3) 
effective chiral lagrangian and hence only act on s-wave partial waves. They 
are able to fit the total cross section data for various strange particle pro- 
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duction reaction channels without introducing resonance states. It remains 
to be seen whether these models can be further improved to account for 
higher partial waves which are definitely needed to give an accurate de- 
scription of the data even at energies near production thresholds. 



3. Unitary ttttN Model 

All of the models described in section 2 rely on the assumption that the 
ttttN continuum can be expanded in terms of quasi-two-particle channels 
such as 7rA, aN , and pN. These models are of course not satisfactory since 
they do not account for all of the effects due to the coupling with the ttttN 
channel. It is necessary to develop a reaction model which also satisfies 
the irirN unitarity condition exactly, This can be done by extending the 
Hamiltonian Eqs.(l)-(3) to include a vertex interaction IV to account for 
the p — ► 7T7T and a —>■ 7T7t decays and to include possible non-resonant tttt 
interaction w TI . Such a formulation and numerical methods for performing 
unitary calculations of two-pion production cross sections are being pursued 
by Lee, Matsuyama, and Sato (LMS) 22 . Here we only briefly describe this 
unitary ttttN model. 

The coupled-channel equations from LMS can also be cast into the form 
of Eqs.(9)-(16) except that the driving term of Eq.(ll) is replaced by 



bq 

'1 & 



V a .b = V n 



X a . b {E) 



(36) 



11 



with 



X a , b (X) = Z a , h {E) + J2 Z a , c (E)g c (E)X C}b (E) . 



(37) 



The driving term of the above integral equation is 



Z a ,b(E) 



i^j 



< a I r+(i) 



p, 



7T7TN 



E - Ho — Vkn.itN - v V7r + ie 



r y (j)|6>(38) 



Note that i ^ j specifing the sum over N* states in the above equation is 
to avoid the double counting of 7T7riV effect which is already included in the 
dressed propagator defined by Eq.(15). 

We have applied this unitary irirN formulation to investigate irN scat- 
tering and 7^V — > irN in Su channel up to W = 2 GcV. The channels 
included are 7r7V, 77.ZV, 7rA and jN. The needed non- resonant interactions 
are generated from tree-diagrams Figs.l(a)-(d) using the unitary transfor- 
mation method. Two N* states are included in the fits to the irN scattering 
amplitude and the E + amplitude of 7./V — > irN. Our results for E + am- 
plitudes are shown in Fig. 2. We see that we are not able to fit the data at 
W > 1.63 GeV and hence only the extracted 7V*(1535) parameters are reli- 
able. Our results are shown in the Table below and compared with the val- 
ues from Chen et al 16 (DMT) and the quark model prediction of Capstick 23 . 
It is interesting to note that LMS's bare value of the iV*(1535) — > 7^ helic- 
ity amplitude A1/2 is close to the quark model prediction. Both the DMT 
and LMS predict that the meson cloud effect, the differences between the 
dressed values and bare values, is to reduce the bare values to the dressed 
values. This is rather different from the situation in the A region where the 
meson cloud is to enhance the transition strength. The differences between 
DMT and LMS values reflect their significant differences in calculating the 
coupled-channel effects. 
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Figure 3. Calculated wN — > ttttN cross sections in S\\ channel. The partial cross 
sections through intermediate irA (pi-D), pN (rho-N) and crN (Sigma-N) are also shown 
to compare with the coherent sum of these channels (Total) , 

To obtain reliable information for the second Sn resonance at about 
1.6 GeV, we are in the process of including pN and aN channels. The 
importance of these two channels can be examined in a unitary calculation 
of irN — > 7T7r7V cross sections. This is achieved by using the Spline-function 
expansion method which was developed in our previous investigations of 
wNN problem. Our results of the partial cross sections of ttN — ► ttttN in 
Sn channel are shown in Fig. 3. Clearly, pN channel must be included for a 
dynamical interpretation of the second N* and to establish whether there 
exists third or even fourth N* in this channel. Our approach is clearly 
different from the investigation of Chen et al. 16 who include only irN and 
rjN channels and the fits to the data are achieved by including up to four 
N*. 



4. Summary 

We have given a unified derivation of most of the models for electromag- 
netic meson production reactions in the nucleon resonance region. An ex- 
tension of the coupled-channel approach to include irirN channel is briefly 
described and some preliminary results for the TV* (1535) excitation have 
been presented. Our complete calculations will be published elsewhere 22 . 
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